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Abstract 

o . 

£NJ t We investigate the asymptotic behavior of the a large class of reversible chemical 

reaction-diffusion equations with the same diffusion. In particular we prove the optimal 
rate in two cases : when there is no diffusion and in the classical "two-by-two" case. 
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1 Introduction 



The self-ionization of water is the chemical reaction in which two water molecules react 
to produce a hydronium ion H%0 + and a hydroxide ion OH~. It can be written as 
follow 

2H 2 =± H 3 + + OH- . (1) 

■ 

This is a classical example of a reversible reaction-diffusion chemical reaction between 
three species such that the water has to appear twice for the reaction. 

We would like to investigate in this paper how fast reversible chemical reactions (as 
the self-ionization) tend to the equilibrium. To be more realistic a diffusion term will 
be added. Such term models the fact that particles have to move (to diffuse) before 
the reaction. If the reaction is not enough mixed, then the diffusion term will be more 
^ ■ important in the equation. 

The main model studied involves q > 1 species Ai interacting together as follows 

g g 

^oaAi^^piAi, (2) 

i=i i=i 

where on, Pi 6 N, with N denoting the set of nonnegative integers. We assume that for 
any 1 < i < q, ai — ^ which correspond to the case of a reaction without catalyzer, 
as in the chemical equation presented in ([T]). 

First we will introduce the differential equations that describe the evolution of the 
species A% for I < i < q in the relation For some < i < q we will denote by ctj the 
concentration of the species A%. We assume existence of two non- negative real- valued 
rate functions 



which describe the evolution of the concentration ai in the two (reversible) reactions. In 
the first one we lose ai molecules of the specie Ai and in the second one there is a gain 
of (3i molecules of the same species A%. We get the opposite for the reverse reaction. 
Thus we can write 

J t a i = (Pi ~ ^ELi o**-*£?=i ft A W 

- (ft - ai)^ELiftA->ELi^A( a )' ( 3 ) 

where a = (01, • • • , a q ). We assume the kinetics to be of mass action type, which means 
that 

9 
3=1 

where I is a positive constant called i/ie rate constant of the reaction. This model was 
proposed by Waage and Guldberg in 1864. The mass action represents the probability of 
the reaction between all the species (a±, Ai)i<i< q which are assumed to be independent. 
Let us denote by k > the rate constant for the reverse reaction. Of course I and k 
could be different. 

We will presume that the pot of the reaction is not mixed or not enough mixed. 
Then concentrations of species depend on the position in the pot and we have to add 
a diffusion term which depends also on the species. We obtain that the model has the 
following mathematical representation 



d t ai = LiOi + (Pi - ai) j I a" J - k J~[ 

3=1 ' 3=1 



q 

a 3 



where for all i, Li is a diffusion operator. 

We refer to [9] for a general introduction on chemical reaction-diffusion models. 



1.1 Mathematical model 

Throughout the entire paper we will assume that all species diffuse with the same speed, 
i.e. for all % £ {1, • • • , q}, Li = L, where L is general diffusion generator. (While this 
is a restrictive case - as generally this hypothesis is not realistic - assuming identical 
diffusion will allow us to obtain some optimal bounds.) 

In more detail we let Q C M. n (n 1) to be a bounded open and connected set and 
assume that the boundary d£l of f2, is C°°-smooth. Let a diffusion operator L be given 
by 

n n 

Lf(x) = Y, ayCxJ^/Cx) +J2h(x)d i f(x), (4) 

i,j=l i=l 

for any smooth functions /, with ajj and b = (6i)i<i< n in C°°(f2) and the matrix 
a(x) = (&ij(x))ij symmetric and positive for all x £ H. (Generally the choice of 
coefficients ajj and bi may depend on the domain SI.) 

We will assume that there exists a probability measure on S7, fdfi = 1, which is L 
- invariant, i.e. for all functions / 6 L 1 (d/i) in the domain of L 

j LfdfjL = 0, 

The existence of an invariant measure is a quite standard problem and has a solution 
for a large class of generators L. 
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If we denote by Oj(t,x), (t ^ 0, x £ f2), the concentrations of the species Ai at time 
t in the position x, then the following reaction-diffusion system is satisfied: 



a 3 



vt > o, vx e n, a^t, s) = Lai + (ft - oi) [ i Yl a°* - k Yl 

. 3=1 J=l 

Vx G 0, aj(0, x) = a°(x) 

Vx e an, vt > o, ^(t,x) = o. 

ov 



(5) 



where k,l > are rate constants of the reverse reaction. It is assumed that the initial 
conditions satisfy for all 1 < % < q, a® ^ and J a®dfi > 0. 

The last equation in ([5]) represents the Neumann boundary conditions which is quite 
natural in the context of chemical reaction-diffusion. 

Let us denote by D(Q) the set of smooth functions / on Q satisfying the Neumann 
boundary conditions, such that for all x £ dfl, j£(t,x) = 0. 

Let (Pt)t^ be the semi-group associated to L, ( which is a linear operator defined 
for all functions / € X>(f2) ), as follows 

vt > o, Vx e n, f P t (/)(x) = LP t (/)(*) 
Vxen, p (/)(x) = /(x) 

Vxean, vt>o, * w; (x) = o. 

ov 

Since // is an invariant measure, so for all t ^ 0, jPtfdfi = J fdfj,. It could be the 
classical heat equation with boundary condition but also the heat equation with a drift 
part some example will be given in the example [TJ 

One of the main tools of the article is the spectral gap inequality (also called Poincare 
inequality). We assume that the operator L satisfies a spectral gap inequality in T>(U), 
that is : There exists Csg ^ such that for all smooth functions / € T>(£1), 

Vai>(/) := J (f- J fdf?j dfi < -2C SG J fLfdfi, (6) 
which is equivalent to, for all functions / E T>(Q) and t ^ 0: 

Var M (P t /) < e'^War^f) . (7) 

In the general form, when L is given by ((3]), then for all function / with Neumann 
boundary conditions, 

r fLfdfi = J V/(x) • (a(x)V/(x))^(x), 

where a(x) is the diffusion matrix of the generator L given in The assumption 
of the spectral gap inequality is fundamental. First able it proves that the semigroup 
(Pt)tssQ is ergodic, that is for all / £ T>(Q), 



t—>oo 



lim P t / = / fdfi, (8) 



in L 2 (<i/x). Moreover this inequality proves that the rate of convergence is exponential 
in L 2 (dfi), it is given by the inequality ([7]). See for instance chapter 2 of [2] of a review 
on semigroups tools and ergodic properties. 
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Example 1 Here are some basic examples on a bounded domain f2. 

One can consider the generator L given in where 6 = and a ^ Xld in £1 with 
A > in the sense of symmetric matrices. If a = Id then the associated semigroup 
is (Pt)t^o * s the classical heat semigroup with Neumann boundary conditions. The 
invariant measure is the Lebesgue measure and the spectral gap constant Csg, depends 
on the domain £1. 

But we can also consider a generator L = A — V^V with i/j a smooth function on 
the domain fL The generator L is symmetric in L 2 (fi^,) where /iw, = e~^dx/C is a 
probability measure in 0, with the normalization constant C . It satisfies a Spectral gap 
inequality with the dirichlet form — J fLfdfj.jp = J|V f\ 2 dix^ for functions f satisfying 
the Neumann Boundary conditions. 



We will first study the problem of a steady state of the differential system © . Let 
Vi(t,x) = \idi(t,x) for some constant Aj > 0, assuming that a\ is solution of ©, then 
Vi satisfies 

Let us chose constants A, > such that kY\i=i(^i) at = 1 111=1 (^«) > or equivalently 
m=i(A) Qi- ^ = We then obtain that Vi is solution of 

Vi > 0, Vx G f2, d t Vi(t,x) = Lvi(t,x) + ki(/3i - ai)G(vi(t,x),- ■ ■ ,v q (t,x)), 

Vx G 0, Vi(0, x) = v®(x) 
ov ■ 

Vx G dn, Vi > 0, -^(t,x) = 
ov 

(9) 

where G(v\, ■■■ , v q ) = ETUi v^ 3 - Ilj=i v j 3 an d 

• v (10) 



Let us define the set S = < (zj)i<j< g , s.t. ^^Zj£;j(/3j — a,f) = >, where fcj is defined 



in (fTUl) . Then for all (zj)i<j< g G 5, one gets 

g g 
dt^ZiVi = L^zm, 

i=l i=l 

which gives 

^2 ZiVi(t, X) = P t I ^ Z i V i ) 
i=l \i=l / 

In particular, due to the fact that fi is an invariant measure, 



g , g , 

Vi(t,x)d/j,(x) = ^Zi J 
i=i i=i ^ 



The goal of this article is to understand the asymptotic behaviour of the reaction- 
diffusion problem in the spirit of [3J 0], for a general but (diagonal) diffusion. 

We will first define a Steady State of equation @, in the following way standard in 
a chemical reversible reaction. 
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Definition 1.1 A steady state of equation (Ejj with non-negative initial conditions 
i v i)l<i<q is a vector i s i)i<i<q £ (R + ) 9 such that for all (zj)i<i<g € S: 

q q n q q 

^2 ^ = ^2 Zi \ and n ^ = n ^ 

i=l i=l J i=l i=l 

Remark 1 We implicitly assume in the previous definition that a steady state is a 
vector (sj)i<j< 9 G independent of x. In fact we don't know in the general case if 

there exists a steady state depending of the space variable x. But as it is proved in this 
paper, the solution converges to the one defined above. 

Lemma 1.2 Let (i>i)i<i< 3 satisfies equation with initial conditions satisfying v® ^ 
and f v^dfj, > 0. Then there exists a unique steady state (sj)i<j<q of @ such that for 
all i G {1, • • • , g} 3 s« > 0. 

Proof 

< The steady state (sj)i<j< ra has to satisfy: 



q q r 

V(zj)i<i<g £ X^ Si = ^Z Zi 

i=l i=l J 



v^dfj, := M z . 



Let ^4 = {(sj)i<j< g , 5Zf=i z i s i = M z }. Since for all i, ct{ ^ then S is a subspace of 
R 9 of dimension g — 1, then A is a manifold of dimension 1. Then one gets 

A=^J v?dii + tki(Pi-ai)J , tGfij. 

We have to find t € R, such that for all 1 < i < q, J v^d/i + tki{j3i — ai) > and 

i / r \ Pi-en 



y / r \Pi—(Xi 

Yli^j v°dn + tkiifa - ai )j = i. 



1=1 

A simple computation gives that the function 



i=i ^ 



is defined on the set [a, 6) with 



a = inf |t, s.t. Vi, ^ w°<i// + tki(/3i — ai) ^ o| 
fe = sup|i, s.t. Wi, J vfd/j, + tk{((3i — ai) ^ Oj . 



For all i G [a, 6), one gets 



2.^ r„p, 



then is increasing and satisfies y(a) = and 93(6) = +00. Thus there exists a unique 
t G (a, 6) such that (p(t) = 1, which ends the proof. l> 
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Due to the fact that all species are moving according to the same diffusion, one can 
reduce the problem as follow: for all 1 < i,j < q, one gets 



and thus 



k(Pi -OLi) ( o ki(fii - Oi) 

ki{fii-ai) 

-vj + Cij{t,x). 



(11) 



Let us fix 1 < i < q. Then the study of ([9]) is equivalent to the study of the following 
PDE with boundary conditions 

d t Vi = Lv-i + ki(Pi - oti) x 

n (H^- + - n (H^. + ^ 

:= Lvi + ki(fii - <Xi)F(t,x,Vi). 

(12) 

The existence problem in the general case when the operators Lj depending on i are 
all different is a difficult problem. Some approaches one can find e.g. in [8j 15} [T2l [TH| [TH] 
and also by discretization in [TT1 [TU] . 

Our contribution to the domain of reaction-diffusion equation is to prove, using 
the spectral gap inequality (or Poincare inequality) and Markov semigroup tools, that 
the solution of reaction-diffusion equation ([9]) converges to the unique steady state 
associated to the initial condition. The goal of this paper is to explain how fast the 
solution converges to the steady state and to debate on the speed of convergence. In 
particular if the result obtained is far or close to the optimal one. 

The main idea is to investigate a simple proof to explain the asymptotic behavior. 
Therefore, for simplicity, we will not focus on existence theorems, the optimal conditions 
of initial conditions and do not study the special case of a unbounded space f2. Almost 
all results given here can be generalize to the entire space M. n if a regular solution of 
the problem is given. 

We will consider different cases. In Section[2]we will study the case without diffusion, 
that means that the concentration of different species do not depend on x G fi. We 
will prove that there exists a solution converging to the unique steady state with an 
exponential and optimal rate of convergence. 

Then in Section [3] we study the classical "two- by- two" case 

A + B^C + V, 

treated in [H [6] by entropy methods. This case is interesting because we obtain the 
optimal rate of convergence and it gives tools to understand a more general case. 
We finish in Section U] with the general case 

q <3 

i=l i=l 

We prove that under the existence of a non-negative solution and under the assumption 
that Oii/3i = 0, we get that solution converges with an exponential rate to the steady 
state. In this case the rate may be not optimal. 
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2 Case without diffusion 



Assume that concentrations of species do not depend on x G Vt. It is the case the pot 
used for the chemical reaction is mixed constantly so that its contents remain spatially 
homogeneous. This case is important because we can solve it explicitly and it gives 
tools to study the general case in the later sections. 

We believe that results of this section are not new, we give here the proof to keep 
the article self contained and also because we did not find appropriate reference. 

The chemical reaction, without diffusion is given by the following system, for all 
ie {l,--- ,q}, 

±vt = kiifii - cm) I fl vf - ft vf | , (13) 
\ i i o i / 



where ki is defined in (jlOp . with initial conditions fj(0) > for 1 < i < q. Equivalently, 
using the same method as for (1121) . for some 1 < i < q fixed, 

d 

dt Vl = 

- ai) (n + c ») - n (M^h + c * 



:= ki(fii - on)F(vi), 



where 



Theorem 2.1 Let the initial conditions (vj(0))i<j< q be positive. Then equation (1 13$ 
has a unique solution defined on [0, +oo), which satisfies for all 1 < i < q, 

\vi(t) - Si\ < e K \vi(0) - Si\ exp (-Ct), 

where K is a constant depending on initial conditions, the steady state (si)i<j<q *s 
defined in Lemma EH and 

c = f[ s ^j2 kM ~ at)2 - ( 15 ) 

• 1 -1 S i 

i=l i=l 

Moreover the constant C is the optimal rate of convergence. 
Proof 

<\ Let jo be such that B~ a — a, n > and for all i s.t. Bi — cti > 0, one has t — /fo — s < 

k ({3i-a ) ■ ^ *he se * 0' A — Q « > 0} is empty one can use the negative part. Assume 
for simplicity that jo = 1. Then the reaction equation becomes 

^-vi = ki(Bi - ai)x 
dt 

By the definition of jo one have, Cn ^ if Bi — on > and Cj i < if Bi — on < 0. 
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Getting a positive solution (i>i)i<j<g) °f is equivalent to getting a solution v\ 
of (|16p which is defined by 

ki(Pi-ai) o ki{Pi-ai) o 

kjiPj-aj) kjify-aj) 3 

then it satisfies the following inequality 

Vte[0,+oo), 0<«i(t)<Af, 

where M = ui(0) — fci(/3i — ax) max |p^z^~yi Pi ~ a i < o|. By convention, if the set 

{i; — a.i < 0} is empty, then we have max j ^jjf^p) i ft — <%i < o| = — oo. 
Let us denote by 

F(X) = k l (p l -a l )x 

Lemma 1 1 . 2 1 proves that the polynomial equation F = has only one solution s\ in the 
set (0, M). Let Q be a factor of F, i.e. we have a factorization F(X) = (X — s±)Q(X). 
Then a simple computation yields 



i=l 



j Jj «,^fe)! := _ c<0) 



i=l i=l 

which proves that F'(si) = Q(si) < and then Q(X) < for all X G (0, M). 

Function F is locally a Lipschitz function, thus by Cauchy-Lipschitz's Theorem, 
there exists a unique maximal solution starting at ^i(O) of the equation 

Vt€[0,r), ^^=F( Ul (t)) = (^(t)- Sl )Q(t; 1 (t)), 

for some T > 0. 

One has Q(X) < for all X G (0, M), which implies that if ui(0) > s±, then 
is non-increasing and moreover vi(t) ^ s\ for all i G [0, T[, while if i>i(0) < si, then 
v\ is non-decreasing and moreover i>i(i) < si for all i G [0, T[, and if «i(0) = s\, then 
= 8l for all i G [0, T[. Then one gets that T = +oo and for all t ^ 0, vi(t) G (0, M). 

Using the identity 

1 _ IjQjsx) R{X) 
F{X) X-8i QiX)' 
with R(X) = (Q(si) - Q(X))/(Q( Sl )(X - si)) we get for all t ^ 

- st\ = \ Vl (0) - si\ exp [Q( Sl )t + y a )- Q( s i) dQ ) . (18) 

This identity proves that v\ goes to si and moreover if 

Jv^o) (a-si)Q(a) 

then for all i ^ 0, 

- si| < el K l|-yi(0) - si\ exp (-Ct). 

Links between V{ and ui gives the last inequalities for any Vi. Moreover, equality (fTKj) 
proves that the constant C is the optimal rate. > 
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3 The "two-by-two" case 



The goal of this section is to investigate the asymptotic behaviour of a chemical reaction 
of particular type 

A + B^± C + V. 

We will assume that all species are moving according to the same generator and the 
speed of the two reactions are the same, for instance equal to 1. 

This case was treated with a general diffusion in [U [6] by using entropy method, but 
the optimal rate of convergence was not obtained. 

Let us denote by a, b, c and d concentrations of A, B, C and T> in the domain Q. In 
this case the functions a, b, c and d are solutions of the following system on J7, 

dta = La — ab + cd 

d t b = Lb-ab + cd , , 

dtc = Lc + ab — cd 
dtd = Ld + ab — cd 

with non-negative initial conditions do, bo, Co and do hi f (^)> such that J a^d^i > and 
the same for bo, Co and do- We recall that the boundary conditions are included in the 
definition of the domain D(JT). 

If a, 6, c, d are solutions, then a + c = Pt(ao + Co), a + d = Pt(ao + do) an d a — b = 
Pt(do — bo) which gives that the function a is solution of the linear equation 

d t a = La - aD t + C t , (20) 

where C t = Pt(a + c )P t (ao + d ) and D t = P t (a + b + c + do). 

Since the solution is given by a linear solution there exists a classical non-negative 
solution of the problem. One can see for example [H [TJ] for a proof . Some remarks on 
the existence are given in Section HI 

We also obtain the estimation useful for the asymptotic behaviour. 

Lemma 3.1 Solution a, b, c and d of (j 19b satisfy for all t ^ and x E 

< a(t,x) < min{P t (a + c ),Pt(a + do)} 
0<b(t,x) < min{P t (6o + c ),Pt(6o + do)} m) 
< c(t,x) < min{P t (co + a ),P t (c + b )} [ ' 

{ 0<d(t,x) <min{P t (d + ao),Pt(do + M} 

Theorem 3.2 Assume that the semigroup (P-t)t^o satisfies a spectral gap inequality Q 
with respect to the invariant probability measure [i. Let ao> bo, cq and do be non-negative 
initial conditions satisfying ao,bo,co,do G L 4 (d/j,). 

i 

We set M a+b+c+d = J(a + b + c + d )d£t, M 4 = (j (a + b + c + do) 4 d^ 2 and 
Csg denotes the constant in the spectral gap inequality Q. 

Let s a is the steady state associated to the initial condition ao, i/M a+ b +c+( i / 8C SG ' 
then the solution a of ()20p satisfies, for all t ^ 0, 



^ J (a - s a ) 2 dfi < (^J J (a - Safd^i + 



5Ma 



^La+b+c+d 8Csg 
exp f- min | M a+b+c+d , ^T^\ tj , ( 22 ) 
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and if M a+b+c+d = then for all t ^ 0, 



sj j (a - s a ) 2 dn < J (oq - Safd^jL + 5M 4 ^ exp (-M a+b+c+d t). 

The same inequality holds for b, c and d associated to s b , s c and s d - 

If the initial conditions satisfy M a+ b+c+d < W~^> ^ en the rate of the convergence 
is optimal. 

Let us start with a general estimate. 

Lemma 3.3 Assume that the semigroup (Pt)t^o satisfies a spectral gap inequality Q 
with respect to the invariant probability measure fi, then for all functions f £ L 4 (ii) and 
all t ^ 0, 



1 -t 



P t / - J fdfij dfi < 4e 2C sg j fHii. (23) 
Proof 

< Set f = f — J fdfx, then using semigroup properties and the Cauchy-Schwartz 
inequality applied to (Pt)t>o ; one gets, 

Pt/- f fdfi) dfi= [ (Pj)%= [ (PtPjYdfi, 



since PtPt/ = Ptf- Now the Markov semigroup (Pt)t^o is given by a Markov kernel, 

2 



so one has ^Pt/J < Pt (f 2 ), which gives at the end 

J (v t f- j fd^j d»< j (?i(?*ff \ dfi. 



If we set F = Pt (/), then one has 



'P t fj 4 dfi<2 j (Pt(F 2 )- / F 2 d/i) d/i + 2( j F 2 d)i 
which gives by definition of F 



4 2^ 



P t /J dn < 2J \ Pt(F*) - J F'dfij dfi + 2^J (P|(/) -Jfdfi) ,/// 

We apply twice inequality (ff]) to F and to / to obtain 

PtfYdfi < 2e"^Var M (F 2 ) + 2e"^ (Var M (/ 2 )) 2 , 



which implies ([23]). > 
Proof of Theorem \SM 

<3 Is this case, the steady state is the following limit, 



s _ J (go + cp) d/i / (go + dp)dfj, _ Hm Ct 
/ (oq + 6q + do + co)<i/Lt t^+oo D t ' 
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the limit can be seen in L 4 {dn). Let us denote by M a+C = J(ao + co)dfj, and define 
similarly M a+d and M a+b+c+d . One has, 

dt \ J Sa ^ d ^ = f( a ~ s a) d tadfi, 
then by ((20|) and ([Mj) one obtains 



c^2~ y (° ~ Safd/s = JaLadfi- M a+b+c+d J (a - s a ) 2 dfi 
+ J a(a- s a )(M a+b+c+d - D t )d/i + J {a - s a ){C t - M a+c M a+d )dfi. (25) 



Let us consider the last term: 



J (a- s a )(C t - M a+c M a+d )dn = J (a - s a )M a+d (P t (a + c ) - M a+c )d^L 

+ j "(a- s a )P t (a + co)(P t (a + do) - M a+d )dfi. 



Setting ip(t) = y f (a — s a ) d/j,, Cauchy-Schwarz inequality yields 



J (a- s a )(C t - M a+c M a+d )dfi < tp(t) \M a+d ^J J (P t (oo + c ) - M a+C ) 2 dfi 
+<p(t)((J p t(«o + c ) 4 df?j 1 Q(P t (ao + d )-M a+d ) 4 d^ 



First spectral gap inequality gives 

y*(Pt(a + c ) - M a+c ) 2 d/j, < e _z? fe*Var M (ao + c ) < e"^* J (a + c ) 2 ^. 
Since the semigroup (Pt)tj>o is contractive : -4 JPt(ao + bo) 4 dfi < 0, one obtains 

(J P t (a + c ) 4 df?j 1 < Q\a + c ) 4 d^ 1 . 
To finish, Lemma 13.31 gives 

(J (P t (oo + d ) - M a+d ) 4 ^ 7 < V2e'^ (/(«o + 4) 4 ^) 7 , 
which implies for the last term of (I25p : 



< 3<p(i) exp ^- J l^J (a + b + c + d ) d\i 
For the other term one gets 

J a(a - s a )(M a+b+c+d - D t )dfjL 

<<p(t)(J a 4 df?j ' (J {M a+b+c+d -D t ) 4 d^j ' . 
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Using (|21j) . we get 



1/4 ( / r \ 1 / 4 / r \ 1 / 4 ' 

^4, ' 



(a + do) a> 



y a 4 dfi^ < min I ^ (a + c ) 4 d^ ' 
and ([23]) gives 

a (a - s a )(M a+b+c+d - D t )dfj, 

i 

<2ip(t)(^j (a + b + co + do^d/A exp ( _ g^;)' 



Then we obtain 



(p'(t) < -M a+b+c+d (f(t) + 5M 4 exp ^ — 8C7gc? ' 

where M4 = f / ( a o + ^o + Co + do) 4 cfy^ ^ ■ Integration of the last differential inequality 

yields: 

if M a 

+b+c+d 7^ 8Csg ' 



I J (a - s a ) 2 d^ < (JJ 



Oo - s a fdn + 



5M 4 



Ma+b+c+d - g^J 



exp j - min <j M a+b+c+d , }l 



and if M a+b+c+d = then 



(o - Safdjl < 



^ j (a - Safdp + 5M 4 t \ exp (-M a+b+c+d t), 



which finished the proof of (|22p . 

If M a+b+c+d < 8 q sg , then the rate becomes e ~ M a+b+c+dt ^ one can c h ec k that M a+b+c+d 
is equal to the constant C defined in (I15D which is optimal. t> 



Remark 2 In the case of a linear equation of diffusion, dtu = Lu where L is given 
by @, the optimal rate of convergence in L 2 (d^i) is given by the spectral gap constant, 
which is the constant Csg i- n inequality ([6]). This rate is independent of the initial 
condition. For reaction- diffusion equations, Theorem \2.1\ and Theorem \3.2\ prove that 
the optimal rate of convergence strongly depends on the initial conditions, which is 
natural for a chemical reaction, the reaction will converges quickly if the species are 
more concentrated at the beginning. 

Remark 3 The result obtained in Theorem VS ."A can be of course generalizes in the case 
when the generator L is just an operator satisfying a spectral gap inequality on the 
domain considered O. For example one can consider the case of a fractional Laplacian, 
a p-Laplacian,..., and others. (A problem may remain to prove the existence of a non- 
negative solution in some of these cases.) 

Remark 4 We recall also that while one can find in the literature some results on 
long time behaviour of solutions, ( as for example nice bounds in TWj), generally it is 
a challenge to obtain the optimal bounds. Contrary to Theorem \3.S\ the rate obtained 
in ffjjl for a general diffusion always depends on the spectral gap constant. 
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4 Study of the general case 



Let us consider now the general case for q ^ 1, 

<? q 

^ CtiA.% ^ — ^ fti-A-i 



i=l i=l 



Assume now that f2 G R n is bounded domain. 

Let consider non-negative initial conditions ^ on fi, for all 1 < i < q. A weak 
solution of ([""J) on the time interval I is q measurable functions (vi) 1<i<q such that for 
all t el, Vi(-,t) G -Li(O), G(vi(-,t), ■ ■ ■ ,v q (-,t)) G Lx(fl), 



• • • ,v q {-,s))\\ Li ds < +oo 
and moreover for all 1 < i < q, x G M n and t > 0, 

v i (t,x)=P t (v?)+ki(p i -a i ) [ Pt- s {G( Vl (;s),--- ,v q (;s)))ds, (26) 

JO 

which satisfies also for all f«(0, x) = ^(x). 

The result given here is a direct application of Rothe [14|, Theorem 4]. The main 
interest of this result is to see why the solution remains non-negative and is defined on 
[0, oo) which generally can be quite surprising for a fully nonlinear parabolic equation; 
(see also [15j for arguments based on a comparison principle). 

Proposition 4.1 Assume that there exist 1 < io,jo < q such that — aj > and 
Pjo ~ <*30 < °- 

Then, for any non-negative bounded and measurable initial condition { v i)i <i<q > 
there exists a non-negative weak solution of the system @ . 

Proof 

<3 We will give here just a sketch of the proof and refer to [H] and references therein 
to get more informations. 

Let us see how solutions are bounded and non-negative. The idea is to solve a 
different problem : As for equation (fTT|) . for all 1 < i < q we note by Ci t \(t,x) = 

P t ( v i ~ ki[pl-a?) v i) ( x ) and consider the PDE 

8 t v 1 = Lv 1 +F{x,t,v 1 ), (27) 

where 



F(x,t,vi) = ki((3i - ai) 

[u {k^cvr + c ^ x) ) + " n l^r^y vi + c ^ x) ) J ' 

and for (x)+ = max{x,0}. 

By [T3], (|27|) has an optimal and weak solution. Let us see why the solution is 
bounded and non-negative. Let v® be a non-negative bounded initial condition and let 
vi a weak solution of ([2*7]) (with the same definition as for equation Q , replacing G by 
F in $26J)). 
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Let for all 1 < i < 



Vi(t,x) = -^4ir — -\vi(t,x) + Cii{t,x). 
ki(Pi - ax) 

Then (v~i)x<i<q is a solution of 

d t Vi = Lvi + ki(/3i - ai)G(vi, ■ ■ ■ ,v g ), 

where 

i=l i=l 

Let us multiply (|29|) by — (v~i)- := min{uj,0}. After integration, we obtain 



(28) 



(29) 



d 



1 



l d\x 



(vi)-L((vi)-)dfi- ki(/3i - Oi) / (vi)-G(vi, ■ ■ ■ ,v q )dfj,, (30) 



where f(vi)-L((vi)-)dfji, = J ViL((vi)-)dfj, < 0. 
On the set {vi < 0}, we have 



h{Pi ~ Ui)G{vx, ■ 



- ai)G(vi, 



.0,- 



where we put at the position i. Since for all j, (vj)+ #J then it is not difficult to see 
that in all cases — a.i)G{v\, ■ ■ ■ , 0, • • • , v q ) ^ 0. Which gives that 



d 
It 



((vi)-) 2 dfi < 0. 



Since at time t = 0, J ((v^)-.) 2 dfj, = then for all t ^ 0, Vi{t) ^ almost everywhere. 

Assume that fi\ — ot\ > 0. Then, since the solutions are non-negative, we get the 
following global estimate of the solution, 



< vi 



kjXPi - ai) 



v jo + 



a 



3o J 



Pt. V 



_ kjXPx - ai) 



'30 



— II 1 Moo 



k\(fix - ax) 



kj (/3 



30 



a 



30 J 



V 



II 

3o Woo' 



If fix — ctx > does not hold, we use io instead jo to get the same result with jo. The 
last estimate proves that vx is bounded and then the solution is defined on [0,oo). 

The same method as above in (|30p proves that for all < i < q, V{ ^ 0. This 
implies that G(vx,-" iVq) = G(vx,--- ,v q ) and then (v~i)x<i<q is also a non-negative 
weak solution of Q which finished the proof of the existence. > 



Remark 5 This restriction on parameters (c^, /3j) is natural in the context of a chem- 
ical reaction by the the principle of conservation of mass by Lavoisier. 

The following corollary is a direct consequence of (jlip . 

Corollary 4.2 Assume that there exist 1 < io, jo < Q such that j3i — ai > and 

fro ~ "jo < 0- 

Let (vi)x<i<q be a solution of ([9]). Then for all 1 < i < q such that /3j — Qj > one 
gets for all t ^ and x G O : 

< »,(*,*) < MA - <*) min {p, (j^L) " y^)) W. " <* <"} 
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and for all 1 < i < q such that — di < : 

< Vi(t,x) < ki\/3i - ai \ mm I P t (-—^ - — — ^ -Vx), fa - ay > o). 

In particular when initial conditions are bounded, solutions of ([5]) are also bounded with 
an explicit upper bound. 

Theorem 4.3 Assume that the semigroup (Pt)t^o satisfies a spectral gap inequality (JZJ) 
with respect to the invariant probability measure [i. Assume also that for all 1 < i < q, 
ctiPi = 0. 

Let (t; 1<i<(? a non-negative bounded initial condition. We assume furthermore 

that for all 1 < i < q, f v®dfi > 0. 

Let (sj)i<i< g be the steady state given by Lemma \1.2l Then for all 1 < i < q, one 
gets 

J J {vi — Si) 2 d/j, < K exp (— min{a, M}t), (31) 

where a > depends on on, /3, and Csg, a nd M, K > depend on the initial conditions. 
Proof 

< The idea is almost the same as in Theorem 13.21 except that we do not obtain the 
optimal rate. 

Assume that f3\ — a\ > 0, the opposite case could be treated in the similar way. 
Equation (I12D applied for j = 1 reads 



d t vi = Lvi + F(t,x,vi), 

where 

F(t,x,y) = - qi)x 

and functions Ci^(t,x) are defined in (fTT]h By the ergodicity properties of the semi- 
group, equation ([8]), implies that in L 2 (dfi), 



lim Ckftz) = / U ~ k h f: "'V i )^ := Q,x 



Denote by F OQ (y) the limit of F(t,x,y) when £ goes to infinity. Note that Foo does not 
depend on x £ H. Then, one gets since J v\Lv\d[i < 0, 



I (vi - sxfdy, < -a^ j ( r - *, )(./•"(/. •. ) - f x (n ))</// 



+k 1 {p 1 -a 1 )j(v 1 -8 1 )F 00 (v 1 )diM. (32) 



Note that -Foo is equal to the polynomial function i 7 defined in (|17|) where in the 
definition of Cy^ in (|14p . Vi(0) is replaced by the mean value of initial conditions J vfdfj,. 
Let us set 

Mx(t, x) = ki(p! - ax) mm jp t (777^ \ " 7777^ H 0*0, Pj ~ <*j < o\, 

\kx(Px-ax) h,i[i ay)/ 
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and for its limit as t goes to oo 

Mr - h(A - «.) min { / (jy^j " JTy^) *. «-«*<«>}■ 

As it was shown in the proof in Theorem 12, 1[ i^fX) = (X —si)Q(X) with Q(X) < 
for all X £ (0,Mf°). Now since for all 1 < i < g, aj/3j = then si is a simple root of 
the polynomial function F^. It implies that that Q(0) < and Q(M^°) < 0. Then by 
continuity of Q there exist e,r/ > such that Q(X) < — e for all X € [0, Mf° + 77]. 

For the second term in ()32|) . we get 

fci(#i. - a\) \ (vx - sx)F 00 (vx)d/j, < -kx(Px - "i)e / («i - «i) 2 d/i 

J y{i)l<Mf°+?7} 



+fci(/3i - «i) / {vx - s 1 )F 00 (v 1 )dn, 

^{ui>Mf +7?} 

and then for some constant -ff depending on initial conditions 

kx(/3x-ax) J (vx - sx)F 00 (v 1 )dfj, < 

- k x {Px -a x )e j (yx - sxfdp, + K/i{vx > Mf° + 77} . 
Corollary 14.21 implies that Mx(t, •) ^ t>i and then Markov inequality gives 

kx{l3i-ax) j (vx - s 1 )F 00 (vx)dfi < 

- kx(Px ~ ai)e J (vx - si) 2 ^ + JTVar^(Mi(t, •)) • (33) 

Since for q measurable functions gi S L^(^) one has 

1 9 

Var M (rnin{#j, 1 < i < q}) < r Var^ft) > 



i=i 

so the last term of ([33|) gives 

kx(Px-ax) J (vx-sx)F 00 (vx)dfi<-kx(^x-ax)e J (vx - Sl ) 2 dfi + K'e'om 1 , (34) 

where K' is an another constant depending on initial conditions. 

Let us note F(t,x,y) = Yll=x K t,i,xV l and F^y) = Y2=x K °o,iy L where we note 
7 = max{^]- =1 fa, Yli=x a i}- The ^ rs ^ term of (1321) gives 

(ui -sx)(F(t , -,Vx) -Poo 



<Jj(vi- sxfd^J J vf(K t ^ x - K^f 



dji. 



i=l 

Let us consider one of them, with 1 < i < g, one has 



J vf{K t ^ x - K^fdfi <JJ vfdpiJ f (K t ,i, x - Kntfdn. 
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There exist some sets Tjj and Aj and constants Hij,k £ H,j,k 6 N such that 
Then for some constant K.% > 0, pj ^ 2 and «• ^ > 0, 

Since the initial conditions are bounded and Pij t k ^ 2, one gets for some another 
constant K 

J ((C lA (t, -)P- k - J ' V- J <ll> < K J (C itl (t, •) - C^fd^ 

and then spectral gap inequality gives for some K' , 



Thus we have proved that there exits 7« > and R{ > depending on initial conditions 
and Csg such that 

All of these estimates give for some a > depending on c^, and Csg and R > 
depending on initial conditions the following bound 



htPt-ai) / («i - st)(F t (vi) - Foaivi^dfi < W / («i - si) 2 dfj,Re 



-at 



Ii(p(t) = y / (01 — si) 2 dfi, the equation (f52"j) becomes, by the previous equation and 

¥>'(*) < -A* (ft - ai)e^(t) + .Re" * + if'e - ^*, 
which finishes the proof. > 



Remark 6 • One can generalize the last theorem in the following way without as- 
suming that for all i, a%Pi = 0. Let us consider i + and j~ such that 



sup 



j,s.t,Pj-aj>a I Pj — a j J Pi+ ~ a i 

and 



sup 



fv°dv\ fv° r dVL 



3,8.t.Pj-otj<0 I Pj ~ a j J Pj~ ~ a j~ 

Assume only that aj+/3j+ = ctj-Pj- = 0. Then the computation for the species i" 1 
and j~ are the same as in the proof of Theorem \J7l 



On can also generalize in assuming that the initial conditions {v®)\<i< q are in 
L q (dfi) for some q > 1 instead of L co (dfi). In that case the proof will be more 
technical. 
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